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Abstract 

Taking the dominant couplings in the standard model to be the quartic scalar coupling, the Yukawa coupling 
of the top quark, and the SU(3) gauge coupling, we consider their associated running couplings to one loop order. 
Despite the non-linear nature of the differential equations governing these functions, we show that they can be 
solved exactly. The nature of these solutions is discussed and their singularity structure is examined. It is shown 
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that for a sufficiently small Higgs mass, the quartic scalar coupling decreases with increasing energy scale and 
becomes negative, indicative of vacuum instability. This behavior changes for a Higgs mass greater than 168 GeV, 
beyond which this couplant increases with increasing energy scales and becomes singular prior to the ultraviolet 
(UV) pole of the Yukawa coupling. Upper and lower bounds on the Higgs mass corresponding to new physics 
at the TeV scale are obtained and compare favourably with the numerical results of the one-loop and two-loop 
analyses with inclusion of electroweak couplings. 



It is well understood that removing infinities arising in quantum field theory when evaluating perturbative cor- 
rections to classical interactions induces a dependence of the physical couplings on the energy scale fx of the process 
being considered. The most surprising consequence of this dependence is asymptotic freedom, the decrease in the 
y—i • non-Abelian gauge coupling with increasing energy scale. Clearly, the asymptotic behaviour of the running couplings 
\ is of great theoretical and phenomenological interest. 

The evolution of these running couplings is dictated by a set of non-linear ordinary differential equations in 
which the various couplings are inextricably linked. This can change a naive expectation of their behaviour based 
on ignoring this linkage. For example, the quartic scalar self-coupling A of a pure O(N) scalar theory diverges with 
I increasing energy scale, but this need not be true if the coupling is affected by the interaction of this scalar with 
other fields. In this paper, we demonstrate via an exact solution that this quartic scalar coupling A in the Standard 
Qh! Model has an asymptotic behaviour and singularity structure that is strongly affected by the top-quark coupling gt 
and the SU(3) gauge coupling 53, the dominant couplings of the Standard Model at the weak scale. We show this by 
explicitly solving in closed form the one-loop equations that govern how A, gt and 53 evolve with fx. The importance 
of having an explicitly analytic solution to a set of non- linear coupling equations has been emphasized in [2]. In fact, 
the asymptotic behaviour and singularity structure of A becomes contingent upon the boundary conditions to these 
equations, and those boundary conditions in turn depend on the Higgs mass Mh- For Mh < 168 GeV we find that 
A decreases with increasing fi and eventually goes negative; for Mh > 168 GeV, A will be shown to remain positive 
up to its singularity. This provides an independent way of using the renormalization group equation to analyze the 
sensitivity of the Standard Model on Mh that is complementary to numerical approaches [4]. Furthermore, having 
an exact solution for these running couplings will provide a way of determining the effective potential at leading- 
log order in the conformal limit of the Standard Model using the method of characteristics pQ, which in turn has 
cosmological implications [5]. 

The one- loop equations that determine how A, gt and g% depend on li are [B]: 
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y = 6y 2 + 3xy - -x 2 
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where x — -S^, y — 



renormalization-induced mass parameter in the theory^ 
Eq. §5§ can be solved immediately to give 



— and the dot denotes the derivative with respect to t — log(/i). Here [i is the 



z(t) = 



l + |z (*-*o) 



where z = z(t ) is a boundary value for z(t). Together, Eq. ([T]) and Eq. (J3j) give 
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dz 
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which, if x — zw, becomes 



dw 
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dz 14 7 
Eq. ([6]) can be immediately integrated to give 
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where wq = w(to) 



fjf^y- If we define K = z ^ 



then Eq. (7) leads to 



x(t) = 



2z(t) 
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We now must examine how the scalar couplant, y, evolves with t. Dividing Eq. (2) by Eq. (3), we obtain 
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From Eq. (7), we can put 
x p 
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t + q 



where t — z ? , p = — ^ and q = — , so that if u — ^ , then Eq. (9) becomes 

t— — - 12u 2 - (6w + 7)u = -3w 2 
dr 

Eq. (Ilip is a Ricatti equation; the transformation 
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u = ^ 

12v 

brings it to the form 
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v = 



where t = Xq and r = = |. A rescaling of v so that 

« = A"(A + 1) /5 / 
leads to 



d 2 f (2a- 
d\ 2 \ A 
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1 It does not appear to be possible to obtain an exact solution by applying computer algebra programs directly to the differential 
equations for these one loop couplings. 
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provided the conditions 
2 25 16 
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are used to eliminate terms proportional to and n+\yi which arise from Eq. (13). Taking the solutions to Eqs. (16) 
and (17) to be 

25 + y/689 
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Eq. (fl~5)) is of the form of a hypergeomctric differential equation whose independent solutions about A = when 
| — A| < 1, are 



f = F(a,b;c;-X) 



and 



/= (-A) 1_c F(a-c+l,&-c+l;2-c;-A) 
provided 
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161 - 3V65 + 3V689 - 13V265 
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a + b = 2a + 2(3 - 7=1 



i(v^689 -V65) 
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c = 2a - 
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1 + -V689 



(24) 



Solving Eqs. ([22]) and (23]) we obtain a = 4 + ±(V689 - V&>) and b = -3 + ±(V689 - v 7 ^). 

The choice of signs in the roots appearing in Eqs. (18) and (19) ensures that c — a — b = — -i- — 2/3 is positive 

so that F(a, b;c;z) is well defined as z — > 1~ as lim 2 ^ 1 - F(a,b; c; z) = ffe~^fe?^Esj ■ As ^ increases, so that — ^ = 

[l + |(i — to)~\ 7 ^1 — exceeds one in magnitude, we must analytically continue the solutions of Eqs. (20, 21) 
using the standard results for analytically continuing the hypergeometric function F(a, b; c; z) into a domain in which 
\z\ > 1. For Re z < |, one can use the result 



F(o, 6; c; z) = (l- z) a F [a,c-b;c 



1-z 



while for |1 — z\ < 1 we have, 

F(a, b; c; z) = ^ c)F(c ~ ° ~ | F (a, 6; 1 + a + & - c; 1 - z) + 
1 (c — a)l (c — o) 

(1 - ^) C - Q - 6 r(C) n ( U?r C) F (« - a, c - M - a - b + c; 1 - z) 

r(a)r(&) 

These extensions are valid for the values of a,b and c given by Eqs. (22-24). However, for the values of — ^ we 
consider, Eqs. (20-21) are valid solutions to Eq. (15). We also note that according to Eq. (8), 

2z(t) 
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so that x(t) develops a pole when — ^ — ► 1. 
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Together Eqs. (12), (14) and (20) show that 
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F(a-c+l,b-c + l;2-c;-z 



12t°(t + qf f F(o, b; c; -Z) + C (-f ) ' ° F(a - c + 1, b - c + 1; 2 - c; -£) ) 
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where, as noted above, r = z 7 , g and C are constants of integration (with q = — ^) and a,b,c,a,(3 are defined by 
Eqs. (18,19,22-24). The derivatives in Eq. (25) can be computed using the equation 



-^-F(a, 6; c; z) = —F(a + 1, b + 1: c + 1: z) , 

az c 

resulting in the expression 
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C(-T/qy- c F(a - c + 2, b - c + 2; 3c; -£) 
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indicating that the Yukawa and scalar couplings share an UV pole at —r/q = 1. In addition, the scalar couplant can 
also develop poles when the denominator shared by the last two terms is zero. Hence, the analytic solution in the 
scalar couplant [Eq. (27)] additionally points to a non-trivial pole structure which builds upon other exact solutions 
in in the literature [31 13] , thereby increasing our understanding of the behaviour of running couplings in the Standard 
Model, most specifically in the case of a scalar field in the presence of both a gauge and Yukawa coupling. 

With initial values of y being taken to be fixed by tree level relation y(<j>o) = „ 2 H ,-i , xq = x{<t>o) = 0.0253, 
zq = z(0o) = 0.0329 [T] (where Mb denotes the mass of a Standard Model Higgs boson and <f>a — 246.2 GeV denotes 
the vacuum expectation value of the scalar field </>), the integration constant C is obtained from Eq. (25) in terms of 



C = 



4.862287086 + 0.00003014398478Mg 
1.236892814 - 0.00004219469394M H 2 



(28) 



where all quantities are in GeV units. 

With these initial values, the energy-dependence of the solution (27) for the scalar couplant y can be examined!! 
For sufficiently small Higgs masses, the couplant decreases with increasing energy and eventually turns negative, 
indicating a vacuum instability in the theory. When the Higgs mass becomes sufficiently large this behavior is seen 
to dramatically change from being a decreasing function to an increasing function of /i, ultimately approaching a 
UV pole. The emergence of a pole in the scalar coupling prior to the pole in the Yukawa coupling indicates that the 
denominators in the last two terms of (28) become zero prior to —r/q = 1. This indicates limits on the applicability 
of perturbation theory for a sufficiently heavy Higgs mass. This is also consistent with the behaviour of the quartic 
scalar couplant reported in [31 HI [7] . 

The value of Mh representing the boundary between these two scenarios corresponds to the case where the pole 
of the Yukawa coupling [and the second term in (27)] coincides with the shared pole of the last two terms in (27). 
In other words, this bound on M# results in a scalar coupling that remains positive until it becomes singular at the 
same energy scale as the Yukawa coupling. Setting the denominator of the last two terms in (27) to zero at the point 
— r/q = 1 leads to the following expression for C 



C 



r(c)r(i - o)r(i - 6) 

T{c-a)T{c-b)T{2-c) 



(29) 



2 We have verified that our analytic results are in agreement with a numerical solution of Eqs. (l)-(3). 
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corresponding via (28) to M H = 168.67 GevH Although the electroweak couplings have been ignored in the RG 
equations (l)-(3), this approximation does not have a significant impact on the numerical value of the boundary 
value of Mh- If one augments the one- loop RG equations to include the electroweak couplings, their numerical 
solution would lead to a value of Mu — 160.8 GeV above which A remains positive up to its singular value. 

Remarkably, the phenomcnological implications of the one-loop analytic expression for the scalar coupling are 
quite consistent with more detailed one-loop and two-loop numerical results including the electroweak gauge cou- 
plings. For example, if we assume a 1 TeV scale at which the Standard Model breaks down, then the naive vacuum 
stability requirement 2/(1 TeV) > leads via (27) to the bound Mh > 75 GeV. By comparison, augmenting the 
one-loop RG equations with electroweak couplings results in Mh > 72 GeV from the numerical solution. Similarly, 
an upper bound on Mh corresponding to a pole at 1 TeV in the last two terms of (27) results in Mh < 740 GeV, 
while the one-loop numerical solution augmented by electroweak couplings results in Mh < 744 GeV. The upper and 
lower bounds on Mh resulting from our analytic solution are also in good agreement with the full two-loop numerical 
analyses [4]. In particular, our bounds on Mh would not appear out of place in Figure 2 of [9 or Figure 3 of [ID] . 

The exact solution of the RG equations in the one-loop approximation with the dominant electroweak scale 
couplings (gt, A, and 53) thus appears to reproduce the essential features of Higgs mass-bound phenomenology. We 
speculate that the underlying singularity structure manifested by the analytic solution (27) is responsible for this 
concordance. 

Finally, we note that knowing the exact behaviour of these one loop couplants would in principle lead to an 
exact expression for the one loop effective potential in the Standard Model when it is computed using the method of 
characteristics [I]. 

Acknowledgements 

NSERC provided financial support. Roger Macleod had helpful advice. 

References 

[1] V. Elias, R.B. Mann, D.G.C. McKeon and T.G. Steele, Nucl. Phys. B 678 (2004) 147; (Erratum) Nucl. Phys. B 
703 (2004) 413. 



[2 
[3 

[4 

[e; 

[7; 

[9; 

[10 



K. A. Meissner and H. Nicolai,"Conformal symmetry and the standard model," arXiv:hep-th/0612165 



M. Lindner, Z. Phys. C31 (1986) 295; 

B. Grzadkowski and M. Lindner, Phys. Lett. B193 (1987) 71; 

B. Grzadkowski, M. Lindner, S. Theisen, Phys. Lett. B 198 (1987) 64. 

J. Kubo, K. Sibold and W. Zimmermann, Nucl. Phys. B259 (1985) 331; Jonathan Bagger, Savas Dimopoulos 
and Eduard Masso, Phys. Rev. Lett. 55 (1985) T. Hambye and K. Riesselmann, Phys. Rev. D 55 (1997) 7255. 

A. Buchel, F. A. Chishtie, V. Elias, K. Freese, R. B. Mann, D. G. C. McKeon and T. G. Steele, JCAP 0503 
(2005) 003. 

M. Slier, Phys. Rept. 179, 273 (1989). 

Joseph D. Lykken, "The standard model: Alchemy and astrology", arXiv:hep-ph/0609274 . 



G. Altarelli and G. Isidori, Phys. Lett. B337 (1994) 141; 

J.A. Casas, J.R. Espinosa and M. Quiros, Phys. Lett. B382 (1996) 374. 

Marcela Carena and Howard E. Haber Prog. Part. Nucl. Phys. 50 (2003) 63. 

Kurt Riesselmann, "Limitations of a Standard Model Higgs Boson", |arXiv:hep-ph/97T l456 



3 



The value Mh = 168.67 GeV obtained from the solution of Eqs. (28) and (29) has also been verified by numerical exploration of 



Eq. (25). 
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